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ABSTRACT 

We used the Near Infrared Camera on Keck I to obtain K s -band images of 
four candidate high-redshift radio galaxies selected using optical and radio data 
in the NOAO Deep Wide-Field Survey in Bootes. Our targets have 1.4 GHz radio 
flux densities greater than 1 mJy, but are undetected in the optical to > 24 Vega 
mag. Spectral energy distribution fitting suggests that three of these objects are 
at z > 3, with radio luminosities near the FR-I / FR-II break. The other has 
photometric redshift z p hot = 1.2, but may in fact be at higher redshift. 

Two of the four objects exhibit diffuse morphologies in iT s -band, suggesting 
that they are still in the process of forming. 

Subject headings: galaxies: active — galaxies: high-redshift — galaxies: starburst 
— infrared: galaxies 

1. Introduction 

High-redshift radio galaxies (HzRGs) are rare objects in the cosmos, residing at the very 
brightest end of the radio luminosity function. Wide-area surveys have identified powerful 
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radio galaxies out to very early cosmic epochs. However, since flux-limited surveys tend 
to probe increasingly luminous objects with increasing redshift, our knowledge of lower- 
luminosity HzRGs at high redshift is sp arse (TN J0924-2201 ha s a 1.4 GHz flux density 
>5i4oo = 73mJy, despite being at z = 5.19; Ivan Breugel et al.lll999l ). We have sought to rem- 



edy t his situation with deep, multifrequency radio maps of the Bootes field (jde Vries et al. 



20021 ; S. Croft et al. (2008, in preparation); §|2]). This field has substantial multiwavelength 
supporting data. 

Many HzRGs sh ow diffuse rest-frame optical morphologies and other characteristics of 
"protogalaxies" (e. g. iPentericci et al. hoOlh : they appear to form when Ly man-break-galaxy- 
sized clumps merge, and eventually evolve into large elliptical galaxies (Ivan Breugel et al. 



19981 ). Identifying lower-luminosity radio galaxies at early cosmic epochs can help our un- 
derstanding of the nature of active galactic nucleus (AGN) hosts as a function of radio 
luminosity, and can help constrain cosmic evolution of the radio galaxy population. 

By applying a cut in radio spectral index, a (S v oc u a ), we can preferentially select 
higher-redshift objects, since a combination of intrinsic evolutionary effects and an observa- 
tional "^-correction" tend to cause a to become ste eper (i. e. , more negative) with increasing 
redshift (for details see, e.g. iKrolik fc Chenlll991i). Ind eed. 35% of sources with a < —1.3 
and 5i4oo > lOmJy observed by lDe Breuck et al.l (120001 ) were found to have z > 3. 



This paper presents X s -band images of radio sources with steeper than normal radio 
spectra and no detections in deep optical data to > 24 mag (Vega) from the NOAO Deep 
Wide-Field Survey (NDWFS) Bootes field, in order to study the morphologies and properties 
of candidate HzRGs that are fainter than "typical" HzRGs. 



Higdon et al.l (120051 ) performed a complementary study of optically-invisible radio sources 



in the Bootes field. They used the Very Large Array (VLA) at 1.4 GHz to map a 0.5 square 
degree area to a limiting sensitivity (at field center) of ~ 15 ^Jybeam -1 . Thirty six of their 
377 radio sources (10%) were not visible in the NDWFS data, although 90% of these 36 
had flux densities less than 1 mJy; in contrast (§ [3]) we consider sources brighter than 1 mJy, 
and with relatively steep spectral indice s, from our ra dio maps (§ covering an order of 
magnitude larger area. Like us, however, iHigdon et al.l conclude that most of t heir optically- 



invisible objects are AGNs at relatively high redshifts. The region mapped by IHigdon et al. 



has very little spatial overlap with our 325 MHz data, and none of their optically-invisible 
radio sources is present in our radio catalog. 



Throughout this paper, we use Vega magnitudes, J2 000 coordinates, an d assume an 
= 0.24, Q A = 0.76, H = 73 km s" 1 Mpc" 1 cosmology JSpergel et alJl2007h . 
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Radio data 



Observations of a 4.9 square degree region of the Bootes field were made with the VLA 
A-array at 325 MHz during four runs from 2003 June-August. The CLEAN beam had a 
FWHM of 5" and the limiting sensitivity was a ~ 150 /xJy. The data analysis and reduction 
will be discussed in detail in an upcoming paper (S. Croft et al. (2008, in preparation)). A 
catalog of radio sources (as opposed to components, i. e. , a d ouble-lob ed rad i o gala xy counts 
as one source) was generated using the method described by lde Vries et al.l (120061 ). 



We matched the 325 MHz catalog (S. Croft et al. (200 8, in preparation)) to the 1.4 GHz 
Westerbork Synthesis Radio Telescope (WSRT) dataset of Ide Vries et al.l (120021 ). which has 
a ~ 28 ^uJy, an ellipsoidal CLEAN beam of 13"x27", and an area of 6.68 square degrees. 
The typical positional accuracy at 1.4 GHz for sources in the high signal-to- noise ratio (S/N) 
regime (which is the case for our sources since their fluxes are > 1 mJy) is 0'.'44 (la). Six 
hundred and fifty three sources (in the 4.9 square degree where the two radio datasets 
overlap) were found to be in common and hence had measured radio spectral indices, a. The 
properties of this matched catalog will be discussed by S. Croft et al. (2008, in preparation). 



3. Sample selection 



To se lect HzRG candidates , we combined deep radio data 
1999h . a deep (B w 



NDWFS 
diameter - 



annuzi &: Dev 



Brodwin et al. 



25.3, R 



[2]) with optical data from the 
24.1, I = 23.6, K ~ 19 mag; 3a, 5" 



120061 ) optical / near-IR survey in Bootes. The magnitude limits for 
point sources are deeper, but HzRGs are likely to be resolved in the ~ 1" seeing of NDWFS, 
so we used 4" diameter optical apertures to determine whether or not a source is identified. 
We used the lBrown et al.l (120071 ) /-band selected catalog, which has improved photometry for 
faint objects compared to the public NDWFS DR3 Source Extractor (SExtractor) catalog. 
In computing photometric redshifts, we measured photometry from the NDWFS optical data 
in 4"diameter apertures at each of the radio positions (§ [6]). The NDWFS Bootes field is 
also covered by deep infrared, X-ray, and other observations (§EJ), although these data were 
not available to us at the time of our Near Infrared Camera (NIRC) observations and were 
not used as part of the selection. 



As noted by lBlumenthal fc Mileyl (119791 ). iTielens et al.l (119791 ). and others, radio sources 
with spectral indices that are steeper than normal (average ~ —0.7) have a higher prob- 
ability of being unidentified in optical surveys, a result that holds true for the combination 
of our radio data with the NDWFS i?-band (Fig. [TJ see S. Croft et al. (2008, in prepara- 
tion) for more details). We also see a similar drop in identification fraction for sources with 
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inverte d (a\i®° > 0.2) spect ra. At S1400 > 1 m Jy, the radio source counts are dominated by 



AGNs ( IHopkins et al.ll2000l ). Pure starbursts with flux densities > 1 mJy must either be at 
relatively low redshift, or have an extremely high star-formation rate, and either ought to be 
easily seen in our deep optical data. One way to preferentially target HzRGs, therefore, is 
to look for radio sources brighter than 1 mJy that are optically faint or invisible, although of 
course even in these objects there may be some contribution to the radio luminosity from a 
starburst component. Deep K-band imaging can pro vide information on the morphologies of 



the HzRGs, and the well- known K — z relation (e. g. , iLilly fc Longairlll984l ; |Pe Breuck et al. 



20021 ; IWillott et al.ll2003l ) can be used to estimate their redshifts. 



Four hundred and thirty of the 434 radio sources with S1400 > 1 mJy in the matched 
325 MHz / 1.4 GHz sample are within the NDWFS imaging coverage. Seventy two of these 
430 have no identifications in the NDWFS optical / IR data (which we define as no detections 
within 4" of the 1.4 GHz radio position in B^RIK to 3cr in a 4" diameter aperture). Some 
of these 72 objects have confusing or extended radio contours, are near bright stars, or ap- 
pear to be associated with extended objects (e. g. nearby galaxies larger than 8" in diameter 
where the radio position is offset from the optical centroid). We inspected finding charts of 
the NDWFS + 325 MHz radio data, and selected a subsample of 14 sources that were un- 
ambiguously unidentified. These 14 sources have radio morphologies that appear unresolved 
at 325 MHz, with no nearby extended radio emission, and are not near bright objects. Only 
one of the 14 has a steep spectrum, with a = —1.48; four have —0.98 < a < —0.87, seven 
have —0.80 < a < —0.59, and two have a > 0.1. The five radio sources with the steepest 
spectra (Fig. [2]) were followed up with infrared imaging (§H]). Here we refer to these sources 
as A - E; see Table [1] for positions and names. 

Despite the efficiency of the method used for combining multiple radio source compo- 
nents in order to eliminate radio lobes from the catalog (§ H]), inevitably a few lobes may 
slip through. One such object turns out to be Source D; it was undetected in our NIRC 
observations, and later turned out to be undetected in any of the IRAC bands. When we 
overplotted our VLA data, as well as data from the NRAO VLA Sky Survey (NVSS), on a 
wide-field Digitized Sky Survey (DSS) image, we found that this source is in fact a lobe of a 
large (4f5 diameter) FR-II that has a faint i?-band identification in the DSS. Additionally, 
after the observations were taken, we replotted the radio data for Source B on the NDWFS 
images, including a lower contour than in our original plots, and noticed that the radio emis- 
sion appeared extended (Fig. [2]), and so this object ought not to have been included in our 
original sample of 14 objects with unresolved radio emission and no NDWFS matches. Nev- 
ertheless, the radio emission appears to be associated with an IR-bright source which meets 
our criteria for optical faintness. Source D (the large FR-II) will not be discussed further 
here; the remaining four fields resulted in the detection of candidate HzRGs, as discussed in 
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Fig. 1. — Fraction of radio sources with S^oo > 1 mJy that are identified in the NDWFS 
i?-band data, as a function of radio spectral index, a. Radio sources with inverted (a > 0.2) 
and steep spectral indices (a < —1) have lower identification fractions than sources with 
"normal" spectral indices, suggesting that they are more distant and / or less luminous. 
This result will be discussed in more detail by S. Croft et al. (2008, in preparation). 



Table 1. Observation log 



Source name 


J2000 ID a 


1.4 GHz RA 


1.4 GHz Dec 


K s R A 


K s Dec 


4.5 fim RA 


4.5 fim Dec 


NIRC exp (s) 


LRIS exp (s) 


Source A 


J142631+341557 


14 26 31.75 


34 15 57.5 


14 26 31.8 


34 15 57 


14 26 31.75 


34 15 58.5 


3840 





Source B 


J142647+345851 


14 26 47.87 


34 58 51.0 


14 26 47.8 


34 58 53 


14 26 47.78 


34 58 53.3 


1920 





Source C 


J142741+342324 


14 27 41.84 


34 23 24.7 


14 27 41.9 


34 23 25 


14 27 41.85 


34 23 26.3 


3840 


3600 


Source D 


J143036+334026 


14 30 36.09 


33 40 26.6 




No ID - radio 


source is a 


lobe 


1920 





Source E b 


J143258+342055 


14 32 58.44 


34 20 55.4 


14 32 58.5 


34 20 56 


14 32 58.44 


34 20 56.0 


1920 


4680 



a J2000 iden tifiers are from the 1.4 GHz catalog of lde Vries et aL | J2002h . All five sources listed correspond to 325 MHz counterparts without ambiguity. The 
Ide Vries et~al] catalog is denoted [DMR2002] by SIMBAD. 

b Detected in XBootes X-ray data (§ 15.11 1 
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Fig. 2.— Postage stamps from NDWFS (B W RIK) and the IRAC Shallow Survey (3.6 
and 4.5 //m) for the five sources in our Keck imaging sample. A 5" radius circle is plotted 
at the 1.4 GHz flux- weighted radio position (jde Vries et al.ll2002l ) in each frame, except K, 
where the 325 MHz radio map is overplotted. The lowest contour is at 0.4 mJy beam -1 ; the 
remaining contours are evenly spaced in flux density up to the peak flux density in each 
map. 
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4. Keck NIRC / LRIS Observations 



Imaging observations were made in i^-band, using the NIRC ( IMatthews fc Soiferlll994j ). 
- a 256 x 256 InSb array with Of 15 pix -1 — on Keck I during two half nights on UT 2004 
April 29 and 30. A 16-position (4 x 4) dither, with dither spacing 3" and 60s per exposure 
(6 x 10 s) was used. Three fields were observed for four dithers (3840 s) and the remaining 
two fields observed for two dithers (1920s), as summarized in Table [U Seeing was Of 5 - Of 6 
and the sky was clear. The data were reduced in the standard manner using DIMSUM0 in 
IRAF. Astrometry was bootstrapped from the NDWFS data; the small offsets between the 
K s and 4.5 jum positions (Table [I]) suggest that our NIRC astrometry is accurate to ~ 1". 
Photometry was performed in a 3f5 diameter aperture; the K s magnitudes are presented in 
Table [21 where we also present the measured NDWFS and ISS photometry (§ED- 

Two of the four HzRG candidates were observed spectroscopically. The slit masks were 



filled w ith secondary targets selected by their IR colors from the AGN "wedge" of IStern et al. 



(120051 ). along with a few HzRG candidates; these will be discussed in f uture papers. T hese 



observations used the Low Resolution Imaging Spectrometer (LRIS; lOke et al.l Il995l ) on 
Keck I, with the D560 dichroic and a slit width of lf3. On the blue side, a 400 line mm -1 
grism, blazed at 3400 A was employed, giving 1.09 A / pixel and spectral resolution 8.1 A. 
On the red side, a 400 line mm -1 grating, blazed at 8500 A was used, giving 1.86 A / pixel 
and spectral resolution 7.3 A. This setup gives spectral coverage of ~ 3150 — 9400 A, with 
some variation from slitlet to slitlet. Source C was observed on UT 2005 June 4 for 3600 s 
(2 x 1800 s) under clear skies, in lf4 seeing. Source E was observed on UT 2005 June 6 for 
4680 s (3 x 1560 s) in Of 8 seeing and with clear skies. The observations are summarized in 
Tabled 

The data were reduced in the standard manner using BOGUS^ in IRAF, and spectra 
extracted in a If 5 wide aperture for many of the secondary targets. The primary targets 
were undetected in both the red and blue side two-dimensional spectra, and for these no 
extractions were performed. 



http : / / iraf .noao . edu/ iraf /f tp/ contrib/dimsumV2 



https : //zwolf kinder . jpl . nasa. gov/$\sim$stern/homepage/bogus .html 



Table 2. Properties of the HzRG candidates 



ID 


„1400 


5i400 (mJy) 


B W a 


R a 


I a 


z b 


J c 


K s d 


[3.6] ° 


[4.5] c 


[5.8] c 


[8.0] c 


Source A 
Source B 
Source C 
Source E 


-1.48 
-0.89 
-0.98 
-0.87 


1.08 ±0.04 
8.93 ±0.04 
3.61 ±0.03 
16.27 ±0.03 


27.60 ± 1.51 
28.25 ± 1.76 
26.70 ± 0.73 
26.37 ±0.76 


25.06 ±0.64 
25.91 ±0.91 
25.32 ± 1.12 
25.73 ± 1.24 


22.50 ±0.29 
27.87 ±3.36 
24.40 ±0.81 
24.48 ± 1.10 


> 23.1 

> 22.9 

> 23.2 


> 20.8 

> 20.6 


20.94 ±0.10 
19.63 ± 0.05 
21.38 ± 0.13 
20.99 ±0.13 


18.12 ±0.14 
18.03 ±0.13 
19.23 ±0.37 
18.45 ±0.19 


17.34 ±0.14 
17.58 ±0.18 
18.51 ±0.40 
17.88 ±0.23 


15.84 ±0.31 
19.10 ±6.05 
17.27 ±1.13 
17.08 ±0.95 


15.98 ±0.57 
15.75 ±0.47 
16.02 ±0.60 
16.42 ±0.85 



Note. — We report measured photometric uncertainties, even when these are large (e. g. 6.05 mag for Source B at 5.8 fJ.va), since these were the values used in photometric redshift 
fitting (§[6). All magnitudes are on the Vega system. 

a NDWFS (§EJ 

b zBo6tes (§[52} 

C FLAMEX (§E3J 

d NIRC (§|U 

°IRAC Shallow Survey (§EU CO 
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Other data 



5.1. X-ray 



The Bootes field is covered by th e Chandra X- r ay ob serva tions of iMurray et all (120051 ). 
The resulting catalogs, published by iKenter et al.l (120051 ) and iBrand et al. reach a 



limiting full-band (0.5 — 7.0 keV) flux of 7.8 x 10 15 erg s 1 cm 2 for on-axis sources with four 
counts. We checked for counterparts to our targets, and found that Source E was detected^] 



(five counts in the full band, corresponding to 1.13 ±0.71 x 10 -15 erg s -1 cm -2 ). This source 
is very likely to be a powerful high redshift a nd / or optically o bscured quasar simply on the 
basis of its large X-ray-to-optical flux ratio (IBrand et al.ll2006l ). 



Intriguingly, an object in the field of Source B is also detected in X-rays with 5 counts 
in the full band (1.17 ± 0.71 x 10 -15 erg s -1 cm -2 ). The X-ray positiorQis close to the peak 
of the 325 MHz radio contours, but offset (by 7'.'0) from the flux-weighted 1.4 GHz position. 
However, it appears that the X-ray and radio sources are two distinct objects; if they were 
the same (or if the radio source were associated with the object nearest the peak of the radio 
contours), the radio source would be extremely core dominated, which is not supported 
by the comparatively steep spectral index. Rather, it appears that the radio source is an 
extended, double-lobed source associated with the J-band dropout at the position given in 
Tabled], and that the X-ray source is associated with the optical object nearest (1'.'2) the X- 
ray position, and clearly visible in the NDWFS J-band (Fig. [2]) and NIRC ii^-band (Fig. [7]) 
data. 

Both of the X-ray sources have too few counts to reliably determine a hardness ratio 
(three counts in the soft band and two in the hard band in both cases). 

The other HzRG candidates are not detected i n XBo otes, although the y lie w ithin the 
coverage area. Since the catalogs of IKenter et al.l (120051 ) and IBrand et al.l (120061 ) contain 
only sources with four or more counts, it is possible that they may be X-ray sources, but 
be insufficiently bright for inclusion in the full catalog. None of our sources is within the 
coverage area of the more sensitiv e (~ 10~ 16 erg s" 1 cm -2 ), single-pointing Chandra X-ray 
observations of I Wang et al.l (120041 ) . 



3 (R.A., decl.) = (14 h 32 m 58.49 s , 34° 20' 54.77") ±1.25" 
4 (R.A., decl.) = (U h 26 m 48.41 s , 34° 58' 47.37") ±1.60" 
5 (R.A., decl.) = (14 h 26 m 48.45 s , 34° 58' 48.47") 
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5.2. Infrared 



The zBootes survey (ICooll 120071 ) provides z-band data for part of the Bootes field. We 
checked for counterparts to our four HzRG candidates, and found that three were undetected 
at the zBootes flux limits (Table [2]), while Source B was outside the zBootes coverage area. 



FLAMEX felston et alihood ) provides J- and -fG-band catalogs for part of the Bootes 



field. The K s -band catalogs are deeper than the K-band catalogs used for our selection, but 
our measured NIRC K s -band fluxes (for the three sources in regions covered by FLAMEX) 
are below the FLAMEX K s -band limit (K s ~ 19.3 mag), so these sources are too faint to 
appear in the i^-band FLAMEX catalog. The J-band FLAMEX limit (J ~ 21 mag) is 
insufficient to put anything but the weakest constraints on J — K s colors, and in any case, 
two of our sources are outside the FLAMEX coverage area. In future, refining our selection 
criteria by requiring non-detections in FLAMEX will further increase our efficiency at finding 
HzRGs, since K > 19 corresponds to z > 2 for sources on the K — z relation (Fig. EJ § [6]). 

W e checked the IRAC ( jFazio et al. I l2004f ) catalogs of the IRAC Shallow Survey (lEisenhardt et al 



20041 ) for infrared counterparts to our HzRG candidates. All were detected at greater than 
2a significance at 3.6 and 4.5 //m, and one also at 5.8 /im. We report the catalog 5" aperture 
magnitudes in Table [21 and use these measurements to assist with spectral energy distribu- 
tion (SED) fitting (§EJ). Table [2] shows measured values, even in cases where the photometric 
uncertainties are large, since real measurements from the data reflect the relative weights 
given to the data points in the \ 2 photometric redshift fitting. All objects are pointlike at 
IRAC resolution (~ 2"). 

We downloaded the post-BCD 24 /im mosaics from the Spitzer archival observations of 
Soifer et al. for the fields surrounding the four HzRG candidates, and performed quick-look 
photometry using the published zero points. Three of our objects were undetected (to a 
flux limit of ~ 0.2 mJy) but Source A (the ultra-steep spectrum source) was detected with 
S 24 /im = 0.43 ± 0.13 mJy. 



6. Discussion 

We have K s and IRAC detections and optical upper limits for our HzRG candidates 
(Table [2]), which we use to make somewhat crude redshift estimates. K s — [3.6], which 
straddles the Balmer and / or 4000 A breaks for the two z p hot > 4 sources from Table [3j 
provides quite tight redshift constraints. The position of our sources on a diagram of [3.6] 
— [4.5] versus K s — [3.6] color (Fig. @J suggests that three are at 3 < z < 5. 
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Fig. 3. — K — z diagram of IWillott et al.l (120031 ) (circles represent radio galaxies with spec- 
troscopic redshifts from the 3CRR, 6CE, 6C*, and 7CRS surveys). Our HzRG candidates are 
plotted as points with error bars, at the redshifts obtained from photometric redshift fitting. 
Vertical error bars represent the measured uncertainties in the K s band magnitudes (where 
we assume K — K s is negligible), and horizontal error bars represent the 68% confidence 
interval in redshift from HyperZ as tabulated in Table 



Table 3. Properties of the HzRG candidates from photometric redshift fitting 



ID 


Zphot 


z (68% confidence) 


z (99% confidence) 


Template 


A v 


M(Brest) 


Rest-frame 1.4GHz L (lO^WHz- 1 ) 


Source A 


4.97 


4.92 - 4.98 


3.94 - 5.69 


Arp220 


0.00 


-24.67 


6.94 


Source B 


3.76 


3.68 - 3.81 


1.77 - 4.08 


QS02 


0.30 


-24.32 


10.5 


Source C 


1.21 


0.00 - 2.76 


0.00 - 6.99 


Mrk231 


0.00 


-19.05 


0.298 


Source E 


4.65 


4.21 - 4.99 


0.00 - 6.99 


Arp220 


0.00 


-24.07 


30.2 



Note. — Precise values of Zp^ot are reported, even when the associated uncertainties arc large, so that the reader may know the exact 
rcdshifts used to compute rest-frame properties of the sources 
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Fig. 4. — [3.6] — [4.5] versus K s — [3.6] diagram for our four HzRG candidates (magenta 
crosses with photometric error bars) , compared to the colors of field sources in Bootes with at 
least a 5<r detection at NDWFS JT-band (where we assume that K — K s is negligible), 3.6 fim 
and 4.5 /im (black po ints). Also plotted are the colors of non-evolving galaxy templates from 
Polletta et al.l (120071 ). redshifted from z = to z = 7 in steps of z = 0.1 (with redshifts 
marked in steps of z = 1.0). The red line is a 2 Gyr old elliptical template, the blue line is 
Arp220, and the green line is a QS02. 
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We also obtained photometric redshifts using the publicly available software HyperZ 
(Version 1.2, iBolzonella et al.ll2000l ). In order to obtain more stringent limits on the fluxes 
of our objects than the 5" aperture limits used in the selection criteria, we remeasured 
photometry from the NDWFS B\yRI images, in a 4" aperture at the position of the radi o 
centroid. Empirical photometric errors were derived using the method of lBrown et al.l ( 20071 ). 
a Monte Carlo technique which places apertures on random regions of sky. This provides 
more accurate estimates of uncertainties (verified by the addition of simulated galaxies to 
the real data) than methods which assume uncorrelated Poisson background noise and no 
source confusion. Although some of the resulting measurements have large uncertainties 
(Table [2]), they help to constrain the resulting photometric redshifts better than using the 
published survey limiting magnitudes as upper limits. In cases where they are large, the 
quoted uncertainties provide a real measure of the variations in the background values due 
to neighboring objects. Since the uncertainties are used essentially as weights in the x 2 
photometric fitting, large uncertainties downweight the effect of bands where the object flux 
is poorly constrained, as opposed to assigning these bands too high a weight by setting the 
fluxes to zero and obtaining the uncertainty from the survey flux limit. However, even when 
we re-run the photometric redshift fitting using the survey limits for those sources with 
uncertainties larger than 1 mag, we find that there is essentially no change in the best-fit 
redshift (the largest change is for Source E, whose best-fitting template changes from Arp220 
to a 2 Gyr old elliptical, and whose redshift changes by just 0.04), although the confidence 
intervals become marginally less well constrained. For the z- and J-bands, where we only 
have access to catalogs and not to the original images (in order to remeasure photometry at 
the radio positions), we set the fluxes to zero and derive the associated uncertainties from 
the published survey limiting magnitudes. We find that whether or not the z- and J-band 
data points are included has essentially no effect on the best-fitting redshifts and confidence 
intervals (due to the large uncertainties in these bands). 

The /-band point for Source A seems rather high (Fig. EJ); inspection of the NDWFS 
/-band postage stamp in Fig. [2] suggests that this is probably due to the unrelated saturated 
star to the so uthwest. In fact, ey en though Source A fulfils our initial selection criterion (no 
source in the iBrown et al.l 120071 catalog closer than 4"), the aperture magnitude measured 
at the radio position (Table [2]) is brighter than the 3a survey flux limits quoted in § El 
strongly suggesting that light from a nearby object is contributing to the measured aperture 
magnitude (the NDWFS /-band point-spread function (PSF) has wings which can cause 
problems for aperture photometry at surprisingly large radii) . If the photometric redshift is 
correct, there could also be some contribution to the flux in this band from Lya at z ~ 5, 
although an / = 22.5 source ought to be visible in the postage stamps (Fig. [2]). Excluding 
the /-band point from the SED fitting for this object has no effect on the best-fitting redshift, 
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and little effect on the confidence intervals. 

Instead of using the templates provided with HyperZ, we used a selection of templates 
(the QS02, Arp220, Mrk23 1, EU2, EU5 and E111 3 templates; the latter three are 2, 5, and 



13 Gyr old ellipticals) from iPolletta et al.l (120071 ). which do a better job of r eproducing IR 



fluxes . The effects of dust extinction were simulated using the prescriptions of ICalzetti et al. 



( 120001 ) . with extinction A v allowed to vary between 0.0 and 2.0 mag (although some of the 
templates already incorporate the effects of some intrinsic dust extinction). B-band absolute 
magnitude, Mb was constrained to be —28.8 < Mb < —19. Because of the faintness of these 
objects in the optical, and the corresponding large uncertainties for the optical data points 
compared to the IR, the range of plausible redshifts is rather large (Table [3]). However, we 
can still estimate the best-fitting redshifts, z p hot- The template SEDs for the best-fitting 
models are shown in Fig. [5J 

Three of the four HzRG candidates have z p hot > 3, and the fourth is rather poorly 
constrained (due to low S/N) to z p hot = 1.2 (Table [3]). We report precise values for z p hot in 
order that the reader may know what values were used for the calculation of the rest-frame 
radio luminosities and other such properties. However, as can be seen from the confidence 
limits in redshift reported in Table [31 we cannot constrain redshifts for these objects very 
accurately using SED fitting, given the large number of bands with non-detections or large 
photometric errors. Source B was not observed spectroscopically, so we cannot say whether 
or not the emission lines seen in the template (Fig. \5§ are really present in its spectrum. In 
the event of galaxy- wide dust extinction these spectral lines would be less prominent, but 
this ought not to affect our determination of z p hot too much. 

We can also use the 4.5 ixm flux densities alone as an alternative to the K — z relation. 



Seymour et all (120071 ) studied a sample of 69 HzRGs with spectroscopic redshift 1 < z < 5.2 
and a wide range of radio luminosities. They showed that the rest-frame i/-band light 
for these objects is well described by models of elliptical galaxies with formation redshift 
Zf = 10 and masses betwee n 10 11 and 10 12 M^. In Fig. [BJ we plot the measured values of 



the 4.5 um flux density fr om lSeymour et all against redshift, along with tracks representing 
Bruzual fe Charlotl (119931 ) models with Zf = 10 and solar metallicity, normalized to 10 11 and 
10 12 Mq. Plotting our HzRG candidates in this figure at their best-fit photometric redshifts, 
these redshifts seem reasonable, although the lower bounds on the redshift (assuming our 
objects have masses > 10 11 M Q ) are still quite weak (z > 2). Source C, with z p hot = 1.21, 
falls below the 10 11 M line, suggesting that it is either unusually low mass, or that the SED 
fitting underestimates the true redshift (certainly possible given the wide confidence interval 
and low S/N for this source). Source C is the only outlier on the K — z plot (Fig. [3]), again 
suggesting it is either very underluminous, or is in fact at higher redshift than suggested by 
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the SED fits. 



The NIRC images (Fig. [7]) provide much higher resolution than IRAC, and we are 
able to see that two of the detected sources have somewhat extended, "fuzzy" morphologies 
(Fig. [8]). Source E is somewhat more compact, but appears to have a faint trail of material 
extending to the northeast. Source B seems to be quite compact and relatively bright. The 
extended sources have if. 9 -ba rid morphologies suggest ive of the optical morphologies seen in 
HzRGs such as B2 0902+34 (Ivan Breugel et al.lll998l ). which hints that these objects may 
be galaxies still in the process of formation (again arguing for a higher redshift for Source C 
than that obtained from SED fitting). The compactness of the observed radio morphologies 
(unresolved in our maps, except for Source B) suggests that we are probably seeing these 
radio galaxies relatively soon after the AGN turn on. The deduced fc-corrected rest-frame 
.B-band magnitudes (Table [3D imply optical luminositi es a factor of a few times fainter than 
the most powerful z > 3 radio galaxies as studied by Ivan Breugel et all (119981 ) and others, 
and the extended emission is a factor of a few smaller (roughly 10 — 30kpc as compared 
to 50 — 100 kpc for the more powerful sources). The extended radio emission of Source B 
(~ 200 kpc if its Zphot is correct) is typical of a classical double-lobed radio galaxy, which 
lends weight to the (independent) selection of a QS02 template in the SED fitting. 



Radio galaxies are among the most massive objects at all redshifts (jSeymour et al.ll2007l ). 
but there is some evidence for a dependenc e of host optical luminosity on radio power (e. g. , 
De Breuck et al.ll2002t iMcLure et all 120041 ) . It is plausible that our objects are less powerful 
in the radio and hence might be expected to be fainter than the K — z relation (Fig. [3D; they 
might, in fact, be still forming the bulk of their stars (as suggested by their morphologies) 
but the AGN has already switched on. Source B, which is more compact in the near-IR, 
and has more extended radio emission (suggesting that the AGN activity is less recent than 
in the other three sources) is slightly brighter than the K — z relation, again as might be 
expected in this scenario. 

From the measured radio fluxes and spectral indices, and the best-fit photometric 
redshifts, we calculate rest-frame 15 1 MHz luminosities (for the z„h. n t > 3 sources) of ~ 
1 x 10 27 W Hz" 1 sr" 1 , which implies Jwillott et al.lll999l : lMcCarthvill993l ) Lya luminosities 
~ 2 x 10 43 ergs _1 . The corresponding Lya fluxes are ~ 6 x 10~ 17 erg s -1 cm -2 , essentially 
independent of redshift (due to similar effects of distance on the observed Lya and radio 
fluxes, and hence the determination of rest-frame luminosities). A typical HzRG Lya line 
(rest-frame FWHM ~ 5 A) would be marginally unresolved in our spectra at z = 0, but 
beyond z ~ 0.6, observed flux densities would start to decrease due to broadening of the 
Lya line in the observed frame. At z ~ 5, the expected Lya flux density is 2 x 10~ 18 erg s -1 
cm -2 A -1 . Since the noise in our spectra has a~4x 10 _19 erg s _1 cm -2 A" 1 , the predicted 
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S/N is rather low (~ 5), so it is perhaps unsurprising that we see no emission lines in our 
spectra, especially given the uncertainties (±1 dex) in the emission line — radio correlation. 
If the line happens to fall in a region of the spectrum where noise is higher (at the position 
of a subtracted telluric emission line, for example), the S/N would be lower still. Dust, if 
present in the source, would also reduce the predicted fluxes, and hence the S/N. Source A 
has the steepest radio spectrum of our candidates, so is most similar to classical HzRGs, and 
remains a good candidate for future spectroscopy. As discussed above, Source B may have 
prominent emission lines, and spectroscopic observations would also be worthwhile. 

The photometric redshifts and "typical" HzRG morphologies support the interpreta- 
tion of these objects as HzRGs (and their deduced radio luminosities and P-band absolute 
magnitudes are AGN-like, as expected). At the fainter radio flux levels we probe, we seem 
to be seeing the more "normal" members of the HzRG population that perhaps make up 
the bulk of the lumino sity function at these redshifts. Indeed, comparing to the results of 



Ledlow fc Owen! (119961 ) , we find that our three z > 3 galaxies (with rest-frame P-band mag- 
nitudes Mr ~ —25.5 mag, as computed from the photometric redshifts in Table [3]) ought to 
lie on the FR-I / FR-II break (~ 10 WHz -1 at these optical absolute magnitudes). 

As a sanity check, we obtain a rough estimate of the space density of these objects by 
calculating the comoving volume covered by the 4.9 square degree combined radio dataset, 
in the redshift range 1.21 - 4.97 (the range of photometric redshifts in Table |3]). In our 
chosen cosmology, this volume is 2.3 x 10 8 Mpc 3 . Since we do not claim to have exhaustively 
studied the radio-loud optically-faint galaxy population in this field within this redshift 
range, the resulting space density for our four HzRG candidates of 2 x 10 -8 Mpc -3 should 
be taken as a lower limit. This compares favorably with the lower limit on radio galaxies of 
si milar radio power (Punn > 10 27 WHz _1 ) but steeper spectral index (a < —1.3) obtained 



by lDe Breuck et al.l (120061 ). 1.2 x 10 9 Mpc 3 at 3 < z < 4. Since only one out of four of our 



HzRG candidates has a spectral index as steep as this (and s ince both estima tes are lower 



limits), the space density is consistent with that measured by iDe Breuck et al 



Deeper optical and / or IR observations (particularly J, H and / or mid-IR) would 
provide better constraints on our photometric redshifts. As discussed above, our objects 
may be only marginally too faint to be detected in our LRIS observations, so it is possible 
that deeper (red-sensitive) 10 m spectroscopy could in future provide spectroscopic redshifts. 
Mid-IR observations would also be of interest in constraining star-formation rates, and the 
presence and temperature of any dust in these objects. 



We thank S. Kulkarni and collaborators for making the LRIS observations of UT 2005 
June 6 (in exchange for images of GRB 0505096 on UT 2005 June 3-4). We thank the MIPS 



-19- 



Bootes team for MIPS data on our sources. We thank M. Polletta and collaborators for mak- 
ing the SWIRE template library publicly available. Data presented herein were obtained at 
the W. M. Keck Observatory, which is operated as a scientific partnership among the Cal- 
ifornia Institute of Technology, the University of California, and the National Aeronautics 
and Space Administration. The Observatory was made possible by the generous financial 
support of the W. M. Keck Foundation. The authors wish to recognize and acknowledge the 
very significant cultural role and reverence that the summit of Mauna Kea has always had 
within the indigenous Hawaiian community. We are most fortunate to have the opportunity 
to conduct observations from this mountain. The work of SC, WvB, WdV and SAS was 
performed under the auspices of the U. S. Department of Energy, National Nuclear Security 
Administration by the University of California, Lawrence Livermore National Laboratory un- 
der contract No. W-7405-Eng-48. SC and WvB acknowledge support for radio galaxy studies 
at UC Merced, including the work reported here, with the Hubble Space Telescope and the 
Spitzer Space Telescope via NASA grants HST #10127, SST #1264353, SST #1265551, SST 
#1279182. The work of DS and PE was carried out at Jet Propulsion Laboratory, California 
Institute of Technology, under a contract with NASA. This work is based in part on ob- 
servations made with the Spitzer Space Telescope, which is operated by the Jet Propulsion 
Laboratory, California Institute of Technology under a contract with NASA. Support for this 
work was provided by NASA through an award issued by JPL / Caltech. This work made 
use of images and data products provided by the NOAO Deep Wide-Field Survey, which is 
supported by the National Optical Astronomy Observatory (NOAO). NOAO is operated by 
AURA, Inc., under a cooperative agreement with the National Science Foundation. Finally, 
we thank the referee, R. Windhorst, for his comments. 

REFERENCES 

Blumenthal, C, & Miley, G. 1979, A&A, 80, 13 

Bolzonella, M., Miralles, J.-M., & Pello, R. 2000, A&A, 363, 476 

Brand, K., et al. 2006, ApJ, 641, 140 

Brodwin, M., et al. 2006, ApJ, 651, 791 

Brown, M. J. I., Dey, A., Jannuzi, B. T., Brand, K., Benson, A. J., Brodwin, M., Croton, 
D. J., & Eisenhardt, P. R. 2007, ApJ, 654, 858 

Bruzual, A. C, & Chariot, S. 1993, ApJ, 405, 538 



-20- 



Calzetti, D., Armus, L., Bohlin, R. C, Kinney, A. L., Koornneef, J., & Storchi-Bergmann, 
T. 2000, ApJ, 533, 682 

Cool, R. J. 2007, ApJS, 169, 21 

De Breuck, C, Klamer, I., Johnston, H., Hunstead, R. W., Bryant, J., Rocca-Volmerange, 
B., & Sadler, E. M. 2006, MNRAS, 366, 58 

De Breuck, C, van Breugel, W., Rottgering, H. J. A., & Miley, G. 2000, A&AS, 143, 303 

De Breuck, C, van Breugel, W., Stanford, S. A., Rottgering, H., Miley, G., & Stern, D. 
2002, AJ, 123, 637 

de Vries, W. H., Becker, R. H., & White, R. L. 2006, A J, 131, 666 

de Vries, W. H., Morganti, R., Rottgering, H. J. A., Vermeulen, R., van Breugel, W., 
Rengelink, R., & Jarvis, M. J. 2002, AJ, 123, 1784 

Eisenhardt, P. R., et al. 2004, ApJS, 154, 48 

Elston, R. J., et al. 2006, ApJ, 639, 816 

Fazio, G. G., et al. 2004, ApJS, 154, 10 

Higdon, J. L., et al. 2005, ApJ, 626, 58 

Hopkins, A., Windhorst, R., Cram, L., & Ekers, R. 2000, Exp. Astron., 10, 419 

Jannuzi, B. T., & Dey, A. 1999, ASP Conf. Ser. 191: Photometric Redshifts and the Detection 
of High Redshift Galaxies, (San Francisco, CA: ASP), 191, 111 

Kenter, A., et al. 2005, ApJS, 161, 9 

Krolik, J. H., & Chen, W. 1991, AJ, 102, 1659 

Ledlow, M. J., & Owen, F. N. 1996, AJ, 112, 9 

Lilly, S. J., & Longair, M. S. 1984, MNRAS, 211, 833 

Matthews, K. & Soifer, B. T. 1994, ASSL Vol. 190: Astronomy with Arrays, The Next 
Generation, 239 

McCarthy, P. J. 1993, ARA&A, 31, 639 

McLure, R. J., Willott, C. J., Jarvis, M. J., Rawlings, S., Hill, G. J., Mitchell, E., Dunlop, 
J. S., & Wold, M. 2004, MNRAS, 351, 347 



Murray, S. S., et al. 2005, ApJS, 161, 1 
Oke, J. B., et al. 1995, PASP, 107, 375 

Pentericci, L., McCarthy, P. J., Rottgering, H. J. A., Miley, G. K., van Breugel, W. J. M., 
& Fosbury, R. 2001, ApJS, 135, 63 

Polletta, M., et al. 2007, |arXiv:astro-ph/0703255| 

Seymour, N., et al. 2007, ApJS, 171, 353 

Spergel, D. N., et al. 2007, ApJS, 170, 377 

Stern, D., et al. 2005, ApJ, 631, 163 

Tielens, A. G. G. M., Miley, G. K., & Willis, A. G. 1979, A&AS, 35, 153 

van Breugel, W., De Breuck, C, Stanford, S. A., Stern, D., Rottgering, H., & Miley, G. 
1999, ApJ, 518, L61 

van Breugel, W. J. M., Stanford, S. A., Spinrad, H., Stern, D., & Graham, J. R. 1998, ApJ, 
502, 614 

Wang, J. X., et al. 2004, AJ, 127, 213 

Willott, C. J., Rawlings, S., Blunder!, K. M., & Lacy, M. 1999, MNRAS, 309, 1017 
Willott, C. J., Rawlings, S., Jarvis, M. J., & Blundell, K. M. 2003, MNRAS, 339, 173 



This preprint was prepared with the AAS IATgX macros v5.2. 



-22 - 





Fig. 5. — SEDs of the HzRG candidates. The best-fitting templates are shown plotted at 
the redshift which gave the lowest xl i n SED fitting. We also note the best-fit redshift in 
each panel; note, however, that the redshifts are not very tightly constrained (Table [3]). 
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Fig. 6. — 4.5 ixm flux density versus redshift for HzRGs with spectroscopic redshifts from 
Seymour et al.l (120071 ) (circles), and our HzRG candidates (crosses, plotted at the best-fit 
photometric redshift with 68% confidence intervals shown). The dashed and dotted lines are 
evolutionary tracks for elliptical galaxies with formation redshift Zf — 10 and masses 10 11 
and 1O 12 M , respectively. 
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Fig. 7. — NIRC Kg-b&nd images of the four HzRG candidates (oriented north up). Overlaid 
are the 325 MHz radio contours, from 0.5 mJy at intervals of lmJy (Source A, Source C), 
from 0.3 mJy at intervals of 0.5 mJy (Source B), and from 0.5 mJy at intervals of 6.0 mJy 
(Source E). To make the K s -band ID more easily visible, the central radio contour is omitted 
in each map. The red circles in the images of Sources B and E mark the positions of the 
XBootes detections (§ 15. ip . 
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Fig. 8. — NIRC images of the four HzRG candidates, centered on the K s positions and 
zoomed to show the details of the radio host -ft^-band morphologies. The exposure times are 
not the same for all four fields (Tabled]). 



